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ABSTRACT 



We report the first simultaneous zJHK spectroscopy on the archetypical Seyfert 2 
Galaxy NGC 1068 covering the wavelength region 0.9 to 2.4/zni. The slit, aligned in the 
NS direction and centred in the optical nucleus, maps a region 300 pc in radius at sub- 
arcsec resolution, with a spectral resolving power of 360 km s""'^. This configuration 
allow us to study the physical properties of the nuclear gas including that of the north 
side of the ionization cone, map the strong excess of continuum emission in the K- 
band and attributed to dust and study the variations, both in flux and profile, in 
the emission lines. Our results show that (1) Mid- to low-ionization emission lines are 
splitted into two components, whose relative strengths vary with the position along 
the slit and seem to be correlated with the jet. (2) The coronal lines are single-peaked 
and are detected only in the central few hundred of parsecs from the nucleus. (3) The 
absorption lines indicate the presence of intermediate age stellar population, which 
might be a signicant contributor to the continuum in the NIR spectra. (4) Through 
some simple photoionization models we find photoionization as the main mechanism 
powering the emitting gas. (5) Calculations using stellar features point to a mass 
concentration inside the 100 - 200 pc of about 10^° Mq. 
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1 INTRODUCTION 

Near-IR spectroscopy has been playing a unique role in our 
understanding of the AGN phenomenon for several reasons. 
First, it includes a wealth of emission lines and stellar ab- 
sorption features not observed in the optical, which are use- 
ful for studying the physical properties of the emission gas 
and to determine the age and metallicity of the nuclear stel- 
lar population as, for instance, the strong CO bandheads in 
the H and K bands (e.g., Schreiber 1998; Origlia, Moorwood, 
& Oliva 1993). Second, extinction by dust is attenuated by 
a factor of ten relative to that of the optical, allowing to 
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probe depths unreachable at shorter wavelengths. Third, it 
is a transition region, where the continuum from the central 
source no longer dominates while the thermal contirmum 
produced by dust and stars becomes important. 

In the last years, the interest for observation in this in- 
terval has increased, and thanks to the availability of cross- 
dispersed (XD) spectrographs, it is possible to study the 
whole 0.8-2.4/im region in a single observation, avoiding the 
aperture and seeing effects that usually affects JHK spec- 
troscopy done in long-slit single band observations. 

With the above in mind, here we present the first spa- 
tially resolved XD spectroscopy covering the interval 0.8- 
2.4/im made for NGC1068. This object is one of the nearest 
and probably the most intensely studied Seyfert 2 galaxy. 
Observations in all wavelength bands from radio to hard X- 
rays have formed a uniquely detailed picture of this source. 
NGC 1068 has played a unique role in the development of 
unified scenarios for Seyfert 1 and Seyfert 2 galaxies (An- 
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tonucci & Miller 1985), in the study of molecular gas in the 
nuclear region of Seyferts (e.g., Myers & Scoville 1987; Tac- 
coni et al. 1994), and in elucidating the importance of star 
formation activity coexistent with the active galactic nu- 
cleus (AGN), on both larger (e.g., Telesco & Decher 1988) 
and smaUer (Macchetto et al. 1994; Thatte et al. 1997) 
scales. NGC 1068 also hosts a prominent narrow-line region 
(NLR) that is approximately co-spatial with a linear radio 
source with two lobes (Wilson & Ulvestad 1983). The narrow 
emission line region has been extensively characterized from 
subarcsecond clouds probed by the Hubble Space Telescope 
(HST) (Evans et al. 1991; Macchetto et al. 1994) to the ion- 
ization cone and extended emission-line region (Pogge 1988; 
Unger et al. 1992) extending to radii of at least 30" (l"=72 
pc at the distance of 14.4 Mpc; TuUy 1988). 

From approximately 2" southwest to the nucleus to 4" 
northeast, emission lines exhibit multiple components (Cecil 
et al. 1990; Crenshaw & Kraemer 2000a). Broad lines were 
found to be approximately 2500 - 4000 km s^^ wide, which 
may be linked to those that are found in polarized light and 
believed to be reflected light from the inner broad line re- 
gion (BLR). Narrow lines consist of a pair of red and blue 
components. Optical studies of [O ll] and [N ll] line profiles 
suggest that the separation of these two components varies 
across the conical NLR. Between ~2".5 and 4". 5 northeast 
from the nucleus, UV line emission is redshifted relative to 
the systemic value, a pattern that is interpreted as the ex- 
pansion of the plasma in the radio lobe (Axon et al. 1998). 

The mechanisms powering the gas of the NLR, i.e. pho- 
toionization from the nucleus or shocks produced by jets, 
have long been under debate. Dopita & Sutherland (1996) 
and Bicknell et al. (1998) proposed that the emission in 
the NLR may be entirely caused by shocks. Velocity split- 
ting over 1000 km^^, reported by Axon et al. (1998), in 
the vicinity of some of the bright emission-line knots pro- 
vides evidence that fast shocks exist in the NLR of NGC 
1068. However, more recent HST data (Crenshaw & Krae- 
mer 2000b; Cecil et al. 2002; Mazzalay et al. 2009) found 
that the emission line ratios are consistent with photoion- 
ization instead of shock heating mechanism. On the other 
hand, models accounting for both photoionization from the 
central radiation and shocks were required to explain both 
the continuum and emission-lines observed in the spectra 
of active galaxies (e.g., Contini, Rodriguez- Ardila & Viegas 
2003, Rodriguez- Ardila, Contini & Viegas 2005). 

The scenario in NGC 1068 is thus clearly complex, and 
observational constraints from all wavelengths need to be 
put together if we want to understand the processes taking 
place in this galaxy. Spatially resolved NLR spectroscopy 
in the NIR can contribute to the solution of the NGC 1068 
puzzle. Here, we present the first observations on this source 
covering simultaneously the 0.8—2.4 /xm region at moderate 
spectral and spatial resolution. Our data are thus able to 
map a wavelength interval rich in emission lines with a large 
range of ionization and absorption features from the stellar 
population not only from the nucleus but also from the cir- 
cumnuclear region, poorly studied in the literature. In §2 a 
description of the observations is given. The main charac- 
teristics of the nuclear and extended spectra are presented 
in §3 and §4 , respectively. Photoionization model predic- 
tions for the most intense lines are discussed in §5, whereas 
the calculations of the central mass based on stellar features 



appear in §6. Finally, our concluding remarks are presented 
in S7. 



2 THE OBSERVATIONS 

The spectra were obtained at the NASA 3m Infrared Tele- 
scope Facility (IRTF) in October 30, 2007. The SpeX spec- 
trograph (Rayner et al., 2003) was used in the short cross- 
dispersed mode (SXD, 0.8 - 2.4 /im). The detector employed 
consisted of a 1024x1024 ALADDIN 3 InSb array with a 
spatial scale of 0.15"/pixGl. A 0.8"x 15" slit oriented in the 
north-south direction was used, providing a spectral reso- 
lution of 360 km/s. This value was determined both from 
the arc lamp and the sky line spectra and was found to 
be constant with wavelength along the observed spectra. 
During the observations the seeing was ~1.3". Observa- 
tions were done nodding in an object-sky-sky-object pat- 
tern with integration time of 180 s per frame and total on- 
source integration time of 12 minutes. After the galaxy, the 
AOV star HD18571 was observed as telluric standard and to 
flux calibrate the object. The spectral reduction, extraction 
and wavelength calibration procedures were performed using 
SPEXTOOL, the in-house software developed and provided 
by the SpeX team for the IRTF community (Gushing, Vacca 
& Rayner, 2004). The S/N of the spectra varies from ~43 at 
the nuclear aperture to ~13 at the most external apertures. 
Intermediate regions have S/N of ~20. 

Figure 1 shows the position of the slit superimposed in 
the galaxy contours obtained from Galliano et al. (2003). 
From our data, we found that the light peak of the nuclear 
aperture in the K-band coincides in position with that deter- 
mined by Galliano et al. (2003) for the same band. The gray 
contours show the 6 cm emission (Gallimore et al. 1996) and 
the red dotted contours show the 20 |J.m image (AUoin et al. 
2000). For NGC 1068, 17 extractions were made along the 
spatial direction: one, of radius 0.8", centered at the peak of 
the K-band distribution, and eight more, of 0.2" in radius, 
at each side of it. In addition, in order to study in detail 
the variations of line profiles and continuum emission with 
distance, the nuclear aperture were subdivided into seven 
smaller apertures (called nuc_a to nuc_g in Figure 1) of 0.2" 
of radius each. We are aware that these smaller apertures 
cannot be used for flux estimates because they are oversam- 
pled. The 1-D spectra were then corrected for telluric ab- 
sorption and flux calibrated using Xtellcor (Vacca, Gushing 
& Rayner, 2003), another in-house software developed by 
the IRTF team. Finally, the different orders of the galaxy 
spectrum were merged to form a single 1-D frame. It was 
later corrected for redshift, determined from the average z 
measured from the most conspicuous lines. No Galactic ex- 
tinction correction was applied because it is negligible at 

NIR (Anir = 0.112). 



3 NUCLEAR SPECTRA 

Below we describe in detail the main features, both in 
emission and absorption, detected in the nuclear region of 
NGC 1068, which correspond to the spectrum covering 1".6, 
centered on the maximum of the galaxy profile. The result- 
ing nuclear spectra in the ZJ, H and K ranges are plotted in 
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Table 1. Observed nuclear emission-line fluxes and EWs 



Line 



A 

(|xm) 



(10- 



-13 



Flux 
ergs cm~ 



-1) 



EW 



[S III] 


0.907 


[S III] 


0.953 


[CI] 


0.985 


[S VIII] 


0.991 


Hell 


1.012 


[SII] 


1.032 


He I 


1.083 


Pa7 


1.093 


[PII] 


1.188 


[SIX] 


1.252 


[Fe II] 


1.257 


Pa^ 


1.282 


[SiX] 


1.430 


[Fe II] 


1.643 


Ha 1-0 S(3) 


1.956 


[Si VI] 


1.962 


H2 1-0 S(l) 


2.121 


Br7 


2.165 


[Ca VIII] 


2.322 



7.50 ± 0.77 
19.80 ± 2.03 
0.91 ± 0.11 
0.91 ± 0.11 
1.61 ± 0.17 
2.81 ± 0.36 
17.00 ± 1.76 
0.97 ± 0.12 
0.78 ± 0.14 
0.88 ± 0.12 
0.91 ± 0.13 
2.24 ± 0.29 
1.59 ± 0.18 
0.72 ± 0.09 
0.54 ± 0.07 
3.39 ± 0.37 
0.35 ± 0.05 
0.42 ± 0.06 
1.01 ± 0.23 



0.0085 
0.0242 
0.0024 
0.0015 
0.0023 
0.0035 
0.0327 
0.0026 
0.0010 
0.0021 
0.0023 
0.0033 
0.0021 
0.0008 
0.0003 
0.0030 
0.0002 
0.0002 
0.0005 



± 0.0008 
± 0.0026 
± 0.0003 
± 0.0002 
± 0.0006 
± 0.0008 
± 0.0071 
± 0.0004 
± 0.0002 
± 0.0003 
± 0.0012 
± 0.0007 
± 0.0007 
± 0.0002 
± 0.0001 
± 0.0009 
± 0.0001 
± 0.0001 
± 0.0002 



Figure 1. NGC 1068 contours obtained from Galliano et al. 
(2003) with the position of the IRTF slit superimposed. The gray 
contours show the 6 cm emission (Gallimore et al. 1996) and the 
red dotted contours show the 20 ^im image (AUoin et al. 2000) 



Figures 2, 3 and 4, respectively, where the wavelengths have 
been translated to the observer's rest frame. 



3.1 Emission lines 

An inspection to Figures 2, 3 and 4 allows us to state 
that the most conspicuous features in the nuclear spec- 
trum are [S ill] AAO. 907,0. 953 ^m, HelAl.083 nm, Pa/3, 
[Fell]A1.643 \mi, Paa, and the coronal lines [Six]A1.430 vim, 
and [Si vi]A1.962 \ua. 'Pary and Br7 are also present, but very 
weak. Pa7 is in the wing of the very strong He I line. Paa 
could not be measured because it partially falls in the gap 
between orders. Emission from H2 molecule is clearly de- 
tected, whose most prominent lines are the transitions H2 
1-0 (S2) and H2 1-0 (SI), although only H2 1-0 (SI) could 
be measured. H2 1-0 (S3) is not clear in the nuclear spec- 
tra because it is too close to the [SiVI] line, which is very 
broad. The flux of H2 can be used to derive the mass of hot 
molecular gas, as shown in section 4.3. 

Several other lines of lower intensity are detected in 
the ZJ spectra, namely [CaI]A0.985 |i.m, [SVIII]A0.991 ^m, 
HenAl.012 ^m, [SII]A1.032 ^m (a blend of 4 [SII] transi- 
tions), [PII]A1.188 ^m, [FeII]A1.257 \mi and the coronal line 
[SIXJA1.252 [im. Wherever possible, emission line fluxes for 
all features were measured assuming that the line profiles 
can be represented by a single or a sum of gaussian profiles. 
The continuum underneath each line was fit by a low-order 
polynomium, usually a straight line. The SMART routine 
(Higdon et al. 2004), a software developed by the Infrared 
Spectrograph (IRS) Instrument Team, originally designed to 
perform real-time processing and analysis of IRS data, was 
used for this purpose. 



The results are shown in Table 1. Errors on the fluxes 
were obtained from the values given by SMART, which is 
calculated based on the residual of the flt and the flux cali- 
bration error, which is about 10 % for the infrared. 

The detection of [P 11] with considerable strength 
is particularly interesting. As explained by Mazzalay & 
Rodriguez- Ardila (2007), in the Sun this element is a factor 
~1000 times less abundant than carbon. In the nuclear spec- 
trum of NGC 1068 we flnd a ratio of [P 11] /[C l] equal to 0.8. 
If the P/C abundance is about solar, the [Pll] line should 
not be present, unless carbon, like iron, might be locked in 
grains. This indicates that phosphorus may be overabun- 
dant by a factor of 10-20 compared to its solar abundance, 
relative to hydrogen. 

The observed line widths are between 800 - 1000 km s~^. 
This is in agreement with what is found in the optical region 
(Pelat & AUoin 1980). Although NGC 1068 has been con- 
sidered a prototypical Seyfert 2 galaxy, its forbidden lines 
are much wider than those of typical Seyfert 2s. Moreover, 
some emission-line proflles show multiple components (Pelat 
& AUoin 1980, AUoin et al. 1983, Meaburn & Pedlar 1986) 
also seen in the NIR as we show in section 4.1. 



3.2 Absorption features 

Several absorption features are clearly seen in the H- and K- 
band nuclear spectrum (Figures 3 and 4). In particular, we 
observe the following ones: Mgl 1.50 |J.m, CO 1.58 |i.m, Sil 
1.59 ^m, CO 1.62 \ira, Mgl 1.71 \ira, CO 2.29 ^m, CO 2.35 
\im and CO 2.38 (xm. These features are identified after their 
main contributor, although, depending on the dominant stel- 
lar type, other features may become dominant. 

The CO 1.62 fi.m can be used to obtain a rough esti- 
mate of the contribution of the stellar population flux to 
the H band. This is done by measuring the depth of the 
^^CO(6-3) overtone bandhead at 1.62 [j.m. We found that 
the observed depth is 12% ± 1% of the continuum. This can 
be compared to the ~ 20% that is typically expected for a 
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Figure 2. Flux calibrated spectrum in the ZJ range of the nuclear region of NGC 1068. The aperture width is 1".6. Major emission 
lines detected are identified. 



population of GKM supergiants that dominates the H band 
light for a stellar population older than 10 yr (Schinnerer 
et al. 1998). It can thus be roughly estimated that about 
60% ± 5% of the total flux in the H band continuum is due 
to GMK supergiants. Following Schinnerer et al. (1998), al- 
though these stars dominate the stellar emission in the H 
band, about one-third of the total stellar flux comes from 
stars of other stellar classes. Therefore, 90% ± 8% of the 
nuclear H band continuum could be of stellar origin. This 
is in reasonable agreement with the 74% ± 8% determined 
by Cid Fernandes et al. (2001) for this object, based on op- 
tical spectroscopy. This means that the H-band continuum 
is essentially dominated by the stellar population, with only 
'-^ 10 - 20 % due to the AGN. This calculations show how the 
NIR spectroscopy is a potential tool to unveil hidden star- 
burst components, either in highly obscured objects or in 
those in which the optical continuum emission is dominated 
by the AGN component (as in Seyfert 1 galaxies). 

One way of estimating the age of the stellar population 
in the NIR is using the CO index defined by Ivanov et al. 
(2000), for the absorption feature starting at 2.29 \ua. This 
index is supposed to be insensitive to extinction and to pos- 
sible uncertainties in the continuum shape of the infrared 
spectra and it is defined as 



CO = -2.5/03 



I (F2.282) 



(1) 



where (-F2.282) and {F2.295) are the average fluxes within 
a bandwidth of 0.01 (i.m, centered on the blue continuum and 
on the bandhead respectively. For the nuclear aperture we 
obtained an index of -0.01 ± 0.05. This value is extremely 
low for a Seyfert 2 galaxy, according to the sample of Ivanov 
et al. Notice that in a Seyfert galaxy nucleus, the late- 



type stellar features are substantially diluted by non-stellar 
nuclear emission. NGC 1068 is probably an extreme case, 
since its nuclear spectrum shows a very steep increase in the 
K band continuum flux with wavelength (Figure 4), much 
steeper than what is expected for a non-thermal power-law 
source. The K band excess of emission is probably due to 
thermal radiation from hot dust (Thatte et al. 1997). 

In order to minimize the influence of the shape of the 
continuum, we normalized the spectrum in this region using 
a linear fit in the region between 2.20 and 2.26 |xm, where 
there are no apparent features. After dividing the spectrum 
by this fit the CO index is recalculated. The new value is 
0.01 ± 0.05. Thus, the low value is not caused by the shape 
of the spectrum and suggests that there is no evidence for 
a strong starburst, in agreement with the results of Ivanov 
et al. (2000) for Seyfert galaxies, but contrary to what was 
found in the H band. 

Nevertheless, it is really hard to estimate the true di- 
lution of the bandhead by the continuum. The best way 
would be to do a full stellar population fitting, including a 
non-stellar component. This is out of the scope of this pa- 
per, but is presented in detail in a follow-up paper (Martins 
et al. 2010). 



4 SPECTRA OF THE EXTENDED REGION 

In addition to the nuclear spectra discussed in the preceding 
section, the spatial information contained in the slit allowed 
us to extract the spectra corresponding to the narrow line 
region and to the extended line region of NGC 1068 to the 
limits of ~350 pc north and south of the nucleus. Gas kine- 
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Figure 3. Same as Figure 2, but in tiie H band. 
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Figure 4. Same as Figure 2, but in the K band. Note tiiat Pa is severely affected by telluric absorption. 



matics as well as extinction and the molecular and [Fen] 
lines will be studied in more detail in the following sections. 



4.1 Gas kinematics 

The spectra for all integrated apertures are presented in Fig- 
ures 5-7, where the wavelengths are in the observer's rest 
frame. It is easy to see that the most conspicuous lines show 
a double peaked profile. One component is in the rest frame 
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of the galaxy (primary component, hereafter) while the other 
(second component, hereafter) is either blueshifted towards 
the north of the nucleus or redshifted to the south. The 
relative intensity of the peaks also varies with distance, as 
shown in Figure 8 for [Slll]A0.953 (xm, the strongest line in 
our spectra. It can be seen that [S ill] is detected at distances 
of up to ~350 pc from the nucleus. As explained in Sec. 2, 
for the purpose of comparing the line profiles we divided 
the nuclear spectra (1.6" aperture) into seven smaller aper- 
tures (see Fig. 1). The aperture named "nuc_a", centered 
on the peak of the continuum light distribution, is extracted 
over the central 0.4" window and the remaining six with 
apertures 0.2" wide. Although the integrated flux of these 
smaller apertures are oversampled because the extraction 
window is smaller than the seeing value, this approach al- 
lows us to study line profiles, which clearly show changes 
along the spatial direction at scales down to 20 pc. 

A comparison of the line profiles in velocity space ob- 
served in the most conspicuous emission lines is presented 
in Figure 9. In this figure, the continuum underneath each 
line was adjusted using a linear interpolation and then sub- 
tracted. It can be seen that the coronal line [Svill] shows 
a profile very similar to that of [Slll], and second compo- 
nents can be clearly seen. The molecular H2 lines also show 
scarce evidence of a second component, not so clear since 
the line is weak. They are also narrower than the emission 
features of other species outside the nucleus, which indicates 
that it might be formed in a different region. Indeed, NIR 
observations has shown that the H2 emission is distributed 
mainly along a direction perpendicular to the ionizing cone 
(Galhano & AUoin 2002). The [Fen] lines seem to follow the 
H2 lines more closely. 

We did not identify a broad component with FWHM 
~ 2000 km s^^ in any line, although they have been re- 
ported on FUSE observations (Zheng et al. 2008) . That may 
be because FUSE slit collects emission from a block of re- 
gions spanning several arcseconds perpendicular to the con- 
ical axis and the total line emission from these regions may 
be blended into one broad component. Since we are observ- 
ing much smaller regions with each aperture, we minimize 
this effect. 

Not all emission lines have double peaked profiles. This 
is the case of high ionization lines such as [S vill] and [Si x] as 
well as molecular lines of H2. Moreover, low ionization lines 
like [Ca l] or [P 11] also seem to be devoided of a double-peak 
structure but these are faint lines. The lack of splitting may 
be either by intrinsic reasons or due to low S/N. In Figure 10, 
the relative flux of the most conspicuous emission lines are 
shown as a function of the distance to the center. The fluxes 
were normalized relative to the nucleus to allow the compar- 
ison between the different lines. Negative distances represent 
the northern direction and positive distances the southern 
one. For double peaked lines, the upper panel shows the dis- 
tribution of the primary component while the bottom panel 
that of the secondary component. It can be seen that the 
second component of [S ill] 0.953/im and Hel 1.083/im are 
observed up to 200 pc from the nucleus. Indeed, it seems 
that there is a maximum of the second components between 
100 - 150 pc from the center. As for the primary compo- 
nents, evidence of emission up to 300 pc from the center 
is observed although it becomes too much faint to be mea- 



sured at 3cr level. Note that the [Fen] emission is the most 
extended of all lines analyzed. 

The flux gas distribution shown in Figure 10 and 
mapped by the first time in the literature with data in 
the interval 0.8— 2.4/Lim confirms previous reports in the 
optical and MIR in the sense that the northern zone has 
stronger emission than the southern one (Mazzalay et al. 
2009; Geballe et al. 2009). This can be attributed to the 
strong extinction towards the observer. In this scenario, the 
ionization cone is partially hidden from our view because it 
crosses the host galaxy plane. Note also that high ionization 
lines like [Six] 1.430^im and [Sivi] 1.954/im are detected 
only up to ~200 pc from the nucleus. This is a clear indica- 
tion that the gas is more ionized near the nucleus. 

The observed emission features of the extended region 
were also measured with the software SMART. The results 
are shown in Tables 2 and 3. Errors were estimated the same 
way as in Table 1. 



4.2 Extinction 

The extinction towards the nucleus and circumnuclear re- 
gion was determined from the comparison of the predicted 
and observed emission line ratios assuming the standard ex- 
tinction law of Cardelli et al. (1989) for Rv = 3.1. From 
our observed lines, two line ratios are extinction indicators: 
Pa/?/Pa7 and [Fen] 1.257^m/1.643^m. We added both com- 
ponents of Pa/3 for this purpose. The reason for this is that 
Pa7 has no second components, probably because it is a 
weak line and it was not possible to distinguish them. When 
using only the primary component of Pa/3, many ratios give 
unphysical results. Paa could not be used in this calcula- 
tion because it was severely affected by telluric transmission. 
Moreover, Br7 was only detected in the central aperture. 

For hydrogen we used the intrinsic ratios for the case B 
given by Osterbrock (1989). The pair of [Fell] lines used to 
estimate the extinction share the same upper level, meaning 
that the intensity ratio depends only on the energy differ- 
ences between the lines and their Einstein- A coefficients, 
making this line ratio a useful probe of reddening. The in- 
trinsic value of the [Fell] 1.257/im/1.643/^m ratio was taken 
from Bautista & Pradhan (1998), calculated from Nuss- 
baumer & Storey (1988) transition probabilities. Table 4 
shows the values of extinction found for each aperture. For 
the errors in this table we only considered the errors in the 
line measurements and not in the flux calibration. 

The results point to extinction variations up to a factor 
of 4 within the NLR. However, the E(B-V) derived from the 
hydrogen ratios suffer from two major problems; the Pa/3 
line may be contaminated with the [Fell] 1.278 (im, lead- 
ing to an overestimation of the reddening determined from 
Pa/3/Pa7. Also, the Pa7 line is in the wing of the much 
stronger Hel 1.083 \x.m line. Therefore, problems with the 
deblending of these two lines could lead to larger uncertain- 
ties in this calculation. Note that the uncertainties quoted 
are determined only from the S/N around the region con- 
taining the lines of interest and they do not take into account 
potential deblending issues. 

Differences between the reddening obtained from the 
[Fen] lines and hydrogen lines are expected, since these lines 
are probably not emitted from the same region. Keep in 
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Figure 5. Spectra of all apertures of NGC 1068 in the ZJ range. From the top to the bottom the spectra correspond to the apertures 
from north to south. 
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Figure 6. Same as in Figure 5, but in the H band. 
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Figure 7. Same as in Figure 5 but in the K band. 
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panels show the apertures to the south, where the second peak is 
second peak is blueshifted. The first plot in each column is the very 



mind that the uncertainties in the E(B-V) from the [Fen] 
lines can also be large since these lines are intrinsically weak. 

The analysis of Kraemer et al. (1998) indicates that 
there may be some dust mixed with gas in various amounts 
in this galaxy. Veilleux et al. (1997) obtained a value of E(B- 
V) = 0.34 using the Balmer decrement. Our average result 
for the NLR from both determinations is above this value. 
This might be an indication that the NIR might come from 
dustier, deeper regions of the NLR. 

The extinction obtained from the two line ratios for 
the various apertures are shown in Figure 11. The nuclear 
region seems to be more dusty than the adjacent regions 
but secondary peaks of extinction are also observed, with 



maximums around 150 - 200 pc from each side of the nucleus. 
These positions are also close to the region where the second 
component of the emission lines has a maximum, and may 
be associated to the presence of dense clouds mixed with 
dust in this region. 

4.3 H2 and [Fe II] lines 

Emission from H2 molecules are observed in most apertures 
of NGC 1068, and has been the subject of extensive anal- 
ysis by different authors. Reunanen et al. (2003) suggest 
a star formation origin for these lines as a result of UV 
heating from OB stars. Whether the heating mechanism is 
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Figure 9. Profile of emission lines in the velocity space. Left panels show the apertures to the south, where the second peak is redshifted. 
Right panels show the apertures to the north, where the second peak is blueshifted. The first plot in each column is the very central 
aperture, plotted in both columns for comparison. Lines were only plotted for the apertures where they have measurable fluxes. 
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Table 2. Observed emission-line fluxes of all South non-nuclear apertures 



South - Main components 


Flux (lO-i"* ergs cm^^ g-i) 


Line 


A (|xm) 


ApOl Ap02 Ap03 


Ap04 


Ap05 


Ap06 


Ap07 


Ap08 



[S III] 

[S III] 

[CI] 

[S VIII] 

He II 

[SII] 

He I 

Pa7 

[PII] 

[SIX] 

[Fe II] 

Pa/3 

[SiX] 

[Fe II] 

[Si VI] 

Ha 1-0 S(l) 

Br7 

[Ca VIII] 



0.907 
0.953 
0.985 
0.991 
1.012 
1.032 
1.083 
1.093 
1.188 
1.252 
1.257 
1.282 
1.430 
1.643 
1.962 
2.121 
2.165 
2.322 



5.81 ± 1.10 
14.21 ± 1.80 
1.60 ± 0.26 
0.21 ± 0.07 

1.12 ± 0.37 
3.50 ± 0.77 
11.48 ± 1.56 

1.13 ± 0.19 
1.10 ± 0.19 
0.39 ± 0.08 
1.39 ± 0.28 
2.22 ± 0.41 
1.43 ± 0.27 
0.76 ± 0.59 
3.28 ± 0.41 
0.65 ± 0.13 

0.21 ± 0.03 



1.63 ± 0.27 
6.66 ± 0.85 
0.48 ± 0.09 

1.00 ± 0.15 
1.28 ± 0.28 

3.64 ± 0.48 
0.74 ± 1.45 
0.82 ± 0.12 
0.11 ± 0.02 
1.06 ± 0.15 
1.04 ± 0.14 
1.25 ± 1.05 
1.27 ± 0.17 
1.06 ± 0.17 
0.32 ± 0.04 



1.28 ± 0.25 
3.95 ± 0.71 
0.45 ± 0.11 

0.63 ± 0.16 
0.81 ± 0.16 
2.36 ± 0.32 
0.58 ± 0.13 
0.64 ± 0.10 
0.72 ± 0.10 

1.06 ± 0.15 

0.57 ± 0.07 
0.17 ± 0.03 
0.24 ± 0.03 



1.31 ± 0.22 
3.93 ± 0.70 
0.64 ± 0.13 



1.04 ± 0.19 0.45 ± 0.08 0.61 ± 0.06 0.23 ± 0.03 
1.84 ± 0.33 0.96 ± 0.14 1.90 ± 0.24 1.22 ± 0.15 



0.61 ± 0.12 0.34 ± 0.10 



2.38 ± 0.43 
0.62 ± 0.14 
1.01 ± 0.13 
0.22 ± 0.09 
0.37 ± 0.08 
0.65 ± 0.12 

0.58 ± 0.08 
0.16 ± 0.03 
0.25 ± 0.03 



0.11 ± 0.02 0.22 ± 0.03 
0.19 ± 0.03 
1.43 ± 0.34 0.77 ± 0.11 0.84 ± 0.12 0.58 ± 0.08 
0.07 ± 0.01 0.08 ± 0.01 



0.28 ± 0.07 0.22 ± 0.05 0.36 ± 0.05 
0.39 ± 0.09 0.13 ± 0.03 0.12 ± 0.02 

0.33 ± 0.06 0.15 ± 0.03 

0.23 ± 0.04 0.12 ± 0.01 
0.14 ± 0.02 0.10 ± 0.02 0.12 ± 0.02 0.10 ± 0.01 



South - Second components 


Flux (10-1"* ergs cm^^ g-i) 


Line 


A (|xm) 


ApOl Ap02 Ap03 


Ap04 


Ap05 


Ap06 


Ap07 


Ap08 



[S III] 


0.907 


2.13 ± 0.44 


2.77 ± 0.52 


1.25 ± 0.13 


1.24 ± 0.24 


[S HI] 


0.953 


4.16 ± 0.81 


4.81 ± 0.87 


2.87 ± 0.31 


2.87 ± 0.54 


Hell 


1.012 


0.17 ± 0.31 


- 


- 


- 


He I 


1.083 


4.09 ± 0.79 


4.46 ± 0.82 


2.20 ± 0.56 


2.18 ± 0.43 


Pa/3 


1.282 


0.20 ± 0.04 


0.29 ± 0.06 


0.18 ± 0.09 


0.76 ± 0.16 


[Si VI] 


1.962 


- 


0.76 ± 0.14 


0.35 ± 0.11 


0.39 ± 0.09 



0.76 ± 0.16 0.65 ± 0.14 

0.69 ± 0.14 0.29 ± 0.07 0.25 ± 0.07 



Table 4. Flux ratios and E(B-V) 



06 05 04 03 02 01 



Apertures 



Line 



Ap07 
Ap06 
Ap05 
Ap04 
Ap03 
Ap02 
ApOl 
Nuc 
Ap09 
AplO 
Apll 
Apl2 
Apl3 
Apl4 
Apl5 



Pa/3/Pa7 



E(B-V) [Fell] 1.257/1.643 E(B-V) 



1.49 ± 0.31 
1.84 ± 0.48 

2.27 ± 0.55 
1.31 ± 0.55 
1.83 ± 0.41 
1.98 ± 0.51 
2.31 ± 0.24 
1.96 ± 0.32 
2.05 ± 0.34 

1.50 ± 0.38 
1.17 ± 0.35 
2.33 ± 0.49 
1.11 ± 0.55 



< 0.01 
0.13 ± 0.25 

1.05 ± 0.26 
0.82 ± 0.56 
0.09 ± 0.22 
0.44 ± 0.17 
1.13 ± 0.10 
0.39 ± 0.12 
0.60 ± 0.13 

< 0.03 

1.18 ± 0.18 

< 0.05 



1.54 ± 0.49 
0.85 ± 0.22 
0.64 ± 0.17 
1.26 ± 0.16 
1.12 ± 0.31 



± 

± 

± 

± 

± 

± 

1.15 ± 0.16 
0.98 ± 0.26 
0.65 ± 0.12 



1.83 
1.26 
1.26 
1.20 
1.09 
0.80 



0.51 
0.15 
0.16 
0.23 
0.10 
0.11 



< 0.92 
1.65 ± 0.26 
2.65 ± 0.26 
0.28 ± 0.13 
0.69 ± 0.28 

< 0.10 
0.26 ± 0.12 
0.28 ± 0.13 
0.44 ± 0.20 
0.80 ± 0.10 
1.90 ± 0.14 
0.60 ± 0.14 
1.17 ± 0.71 
2.60 ± 0.19 




Figure 11. E(B-V) derived from emission lines as a function of 
distance from the center. Negative distances represent the north 
.direction and positive distances the south. 



thermal or fluorescence is not clear because neither the 2- 
1 S(l) 2.247/im nor the 1-0 S(0) 2.223pm emission were de- 
tected in the nucleus. They derived an upper limit of 5300 K 
for the vibrational temperature in the inner few hundred of 
parsecs of this object. 



We used their method to estimate the mass of the 
hot molecular gas present in each aperture. Assuming T = 
2000K, a transition probability Ag(l) = 3.47 x lO"'^ s"^ 
(Turner et al. 1977) and the population fraction in the f = 
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Table 3. Observed emission-line fluxes of all North non-nuclear apertures 



North - Main components 






Flux (10-1"* 


ergs cm ^ s ^) 














Line 


A(|j.m) 


Ap09 


AplO 


Apll 


Apl2 


Apl3 


Apl4 


Apl5 


Apl6 


[S III] 


0.907 


17.96 ± 2.06 


11.73 ± 1.18 


6.50 ± 0.96 


2.53 ± 0.33 


1.74 ± 0.24 


1.43 ± 0.19 


1.32 ± 0.18 


0.56 ± 0.09 


[S III] 


0.953 


46.14 ± 6.96 


30.03 ± 3.01 


16.33 ± 2.10 


7.46 ± 0.99 


4.83 ± 0.62 


2.94 ± 0.39 


3.25 ± 0.45 


1.24 ± 0.15 


[CI] 


0.985 


2.38 ± 0.42 


1.85 ± 0.27 


1.39 ± 0.20 


0.88 ± 0.13 


0.53 ± 0.08 


0.55 ± 0.07 


- 


- 


[S VIII] 


0.991 


1.86 ± 0.30 


1.40 ± 0.21 


0.67 ± 0.10 


0.21 ± 0.04 


- 


- 


- 


- 


Hell 


1.012 


4.25 ± 0.70 


2.19 ± 0.30 


1.09 ± 0.16 


0.68 ± 0.10 


0.50 ± 0.08 


0.57 ± 0.08 


- 


- 


[SII] 


1.032 


6.36 ± 1.05 


4.15 ± 0.57 


2.30 ± 0.30 


0.90 ± 0.12 


0.69 ± 0.09 


0.30 ± 0.05 


- 


- 


He I 


1.083 


29.19 ± 4.34 


17.39 ± 1.75 


9.87 ± 1.33 


4.75 ± 0.65 


2.31 ± 0.30 


2.19 ± 0.30 


1.98 ± 0.22 


0.74 ± 0.13 


Pa7 


1.093 


2.46 ± 0.31 


1.45 ± 0.20 


1.10 ± 0.19 


0.72 ± 0.13 


0.30 ± 0.05 


0.38 ± 0.05 


- 


- 


[PII] 


1.188 


2.01 ± 0.22 


1.46 ± 0.19 


1.01 ± 0.13 


1.01 ± 0.13 


0.39 ± 0.06 


- 


- 


- 


[SIX] 


1.252 


2.31 ± 0.28 


1.68 ± 0.23 


0.81 ± 0.12 


0.99 ± 0.16 


0.36 ± 0.06 


0.22 ± 0.04 


0.15 ± 0.03 


- 


[Fc II] 


1.257 


2.40 ± 0.29 


2.06 ± 0.28 


1.95 ± 0.27 


1.14 ± 0.19 


1.15 ± 0.17 


0.77 ± 0.10 


0.51 ± 0.07 


- 


Pa/3 


1.282 


4.81 ± 0.70 


2.44 ± 0.32 


1.02 ± 0.42 


0.85 ± 0.12 


0.52 ± 0.08 


0.27 ± 0.04 


0.24 ± 0.04 


- 


[SiX] 


1.430 


2.41 ± 0.30 


1.06 ± 0.15 


0.76 ± 0.11 


0.45 ± 0.09 


- 


- 


- 


- 


[Fc II] 


1.643 


1.91 ± 0.26 


1.71 ± 0.25 


1.80 ± 0.21 


1.43 ± 0.16 


1.00 ± 0.12 


0.79 ± 0.09 


0.78 ± 0.10 


- 


[Si VI] 


1.962 


5.51 ± 0.69 


2.42 ± 0.36 


0.86 ± 0.12 


0.43 ± 0.10 


0.13 ± 0.03 


- 


- 


- 


Ha 1-0 S(l) 


2.121 


0.76 ± 0.10 


0.53 ± 0.06 


0.32 ± 0.04 


0.21 ± 0.04 


0.13 ± 0.02 


0.08 ± 0.01 


0.15 ± 0.03 


- 


Br7 


2.165 


0.99 ± 0.12 


0.58 ± 0.09 


- 


0.16 ± 0.04 


- 


0.12 ± 0.01 


- 


- 


[Ca VIII] 


2.322 


1.14 ± 0.17 


0.39 ± 0.09 


- 


- 


- 


- 


- 


- 


North - Second components 






Flux (IQ-i"* 


ergs cm ^ s ^) 














Line 


A (|xm) 


Ap09 


AplO 


Apll 


Apl2 


Apl3 


Apl4 


Apl5 


Apl6 


[S III] 


0.907 


2.49 ± 0.33 


1.66 ± 0.22 


1.20 ± 0.16 


2.29 ± 0.33 


1.18 ± 0.18 


0.55 ± 0.08 


0.43 ± 0.07 


_ 


[S III] 


0.953 


5.01 ± 0.70 


3.90 ± 0.52 


3.07 ± 0.45 


3.33 ± 0.44 


1.89 ± 0.28 


1.20 ± 0.18 


0.88 ± 0.13 


- 


[S VIII] 


0.991 


0.47 ± 0.07 


- 


- 


- 


- 


- 


- 


- 


He II 


1.012 


- 


- 


- 


- 


- 


- 


- 


- 


He I 


1.083 


4.23 ± 0.58 


3.64 ± 0.49 


1.98 ± 0.27 


1.99 ± 0.29 


2.07 ± 0.23 


0.58 ± 0.09 


- 


- 


Pa/3 


1.282 


0.58 ± 0.09 


0.53 ± 0.08 


0.62 ± 0.09 


- 


0.19 ± 0.31 


0.15 ± 0.03 


- 


- 


[SiX] 


1.430 


0.34 ± 0.05 


0.25 ± 0.04 


- 


- 


- 


- 


- 


- 


[Fe II] 


1.643 


0.47 ± 0.06 


0.46 ± 0.07 


- 


- 


- 


- 


- 


- 


[Si VI] 


1.962 


2.69 ± 0.35 


- 


- 


- 


- 


- 


- 


- 



1, J=3 level L=i,j=3 = 0.0122 (Scoville et al. 1982), the H2 
mass was obtained from 



Table 5. Mass of the excited molecular hydrogen 



m//„ 



5.0875a;10^^L>^/i_o5(i)10' 



0.4277^22 



(2) 



where D is the distance in Mpc, Ii-os(i) is the observed 
flux of H2 A2.121 \uii in erg cm^'^ s^'^ and A22 is the 2.2-|im 
extinction calculated through the E(B-V) from Table 4. The 
mass mH2 is given in Mq. 

The results for each aperture are given in Table 5. We 
have to point here that it is now well known that the distri- 
bution of H2 (perpendicular to the NLR) is so that wc miss 
a large fraction of it in our slit (Sanchez et al. 2009, Galliano 
et al. 2003) . These values are in general agreement with the 
values obtained by Reunanen et al. (2003) for Seyfert 1 and 
2 galaxies, although the nuclear value we obtained is much 
smaller than what they obtained for NGC 1068. Several rea- 
sons may be invoked to explain the difference. The aperture 
they use to do the analysis is much larger than ours, en- 
compassing all our apertures. As can be seen in Table 4, the 
extinction varies and the higher values are not in the center. 
Reunanen et al. probably got an average of these values. In 
addition, they obtained A22 from the continuum extinction. 



South 


VCIH^ [Mq] 


North 


iiiHa [M0] 


Ap08 


11.8 


Apl6 


- 


Ap07 


13.5 


Apl5 


49.8 


Ap06 


16.8 


Apl4 


15.1 


Ap05 


30.9 


Apl3 


189.5 


Ap04 


82.0 


Apl2 


65.7 


Ap03 


31.4 


Apll 


52.0 


Ap02 


48.7 


AplO 


72.6 


ApOl 


78.1 


Ap09 


98.2 


Nuc 


449.8 







assuming that the difference for the line extinction would 
not be large. This may be true for many Seyfert galaxies, 
but the continuum of the nuclear region of NGC 1068 seems 
to be extremely dusty, as can be seen in Figure 7, by the 
inclination of the spectrum. 

Using the extinction we derived from the emission lines 
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Figure 10. Relative flux of main emission lines as a function of the distance to the center. The fluxes were normalized by the center 
values to allow comparison between lines. Negative distances represent the north direction and positive distances the south. The top 
panel shows the main components, while the bottom panel shows the second component if present. 



(section 4.2) we obtained a much smaller value for the hot 
molecular gas mass. 



4.4 Excitation mechanisms of the H2 and [Fe II] 
lines 

Galliano et al. (2003) analyze the possibility of explaining 
their observations of the MIR emission with diffuse dust in 
the ionization cone. This model is discarded because the 
diffuse dust would shield the central engine as well as the 
BLR from any direction, requiring an ionizing mechanism 
other than photoionization to explain the NLR emission. 
On the other hand, they favor a model of thick dust clouds 
in the cone, with a small covering factor, heated by the cen- 
tral source radiation, without discarding additional heating 
mechanisms by shocks generated by the interaction of radio 
jets and gas clouds. 

In order to explain the observed molecular emission per- 
pendicular to the cone axis, Galliano et al.(2003) built a nu- 
merical model combining a molecular disc and an axisym- 
metrical, no coplanar. X-ray absorber (see their Figure 9), 



reached by X-rays from the central source. This two compo- 
nent model reproduces quite well the main features observed 
in the H2 and CO emission maps. 

More recently, the excitation of near infrared H2 and 
Fe II lines have been discussed by Rodriguez-Ardila et al. 
(2004) based on observations covering the inner 300 pc 
of most Seyfert galaxies and minimizing contamination by 
the host galaxy. They showed that the diagnostic diagram 
H2A2.121/Br7 versus [FeII]A1.257 /Pa/3 proposed by Larkin 
et al. (1998) is an efficient tool to separate objects by the 
level of nuclear activity, i.e., AGN, Starburst or Liner. By 
including more objects, Rodriguez-Ardila, Riffel & Pastor- 
iza (2005) further test the efficiency of the diagram. For 
NGC 1068, we obtained 0.7 < H2A2.121/Br7 < 1.3 and 0.6 
< [FeII]A1.257 /Pa/3 < 2.8. Our resufis are shown in Figure 
12 (filled circles) together with those observed by Rodriguez- 
Ardila 2004 and 2005 (crosses); the dashed lines show the 
limits suggested by these authors. The observational data 
for NGC 1068 corresponding to the nucleus, as well as to 
different locations along the slit, are in the AGN region of 
the diagram. This reinforces the results found previously, 
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Figure 12. Ratios H2 2.121/Br7 versus [Fe II] 1.257 /Pa/3. Filled 
circles represent our data for NGC 1068. Crosses to data obtained 
by Rodriguez- Ardila et al. (2004, 2005). 



i.e., that this diagram is a helpful tool for classifying ob- 
jects according to their degree of excitation, and shows that 
line ratio intensities from slits at different distances from the 
center can be also used. 

In addition, the [FeII]A1.257/[PII]A1.188 line ratio can 
be used as a tracer of the ionization mechanism of the emit- 
ting gas and depends only on the Fe/P abundance. Following 
Oliva et al. (2001), in photoionized gas, like HII regions and 
planetary nebulae, this line ratio is less than 2, while in shock 
dominated gas its value is higher than 20. The high value 
would be due to the evaporation of the iron-based grains in 
the shock front, increasing the Fe abundance relative to P in 
the low-ionization region. Our data show 0.70 ^ [Fe II]/[P 
II] ^ 3.0. Such low values might indicate that, either like 
we mention in section 3.1, the high value of [PII] observed 
here might be due to an overabundance of phosphorus, or 
most of the Fe is locked in grains and the emitting gas is 
photoionized. 



5 PHYSICAL CONDITIONS OF THE GAS: 
NIR CORONAL AND LOW-IONIZATION 
EMISSION-LINES 

In the last few years, photoionization models have been 
largely used to analyze the physical conditions of the NLR 
of NGC 1068. Some of these studies tried to reproduce the 
observed emission-line intensities in a wide wavelength range 
(Alexander et al. 2000, Kraemer & Crenshaw 2000a, Krae- 
mer & Crenshaw 2000b, Spinoglio et al. 2005), while others 
were based on selected lines (Groves et al. 2004, Groves et al. 
2006, Martins, Viegas & Gruenwald 2003, Rodriguez- Ardila 
et al. 2006). Recently, Osaki (2009) focus his analysis on the 
kinematic and excitation structure of the NLR of NGC 1068. 



Although the basic photoionization codes used in the 
literature are very similar, the choice of the input param- 
eters by different authors may change - in particular, the 
choice of the filling factor. Since one-dimensional photoion- 
ization codes are unable to properly deal with an inhomoge- 
ncous gas distribution, the optical depth at each cloud shell 
is multiplied by a chosen filling factor (see, for instance, Os- 
terbrock 1989). The net effect is to increase the size of the 
low ionized zone as well as to decrease the average tempera- 
ture of the cloud, which could mislead the conclusions. Fol- 
lowing Martins et al. (2003), in the present paper we adopt a 
multicloud model which should provide a better simulation 
of the inhomogeneous NLR of NGC 1068. 

Another important issue concerning photoionization 
models is the use of the so-called ionization parameter U 
to characterize the intensity of the ionizing radiation con- 
tinuum. In fact, U is the ratio of the density of the ionizing 
photons reaching the cloud to the gas density. It is a con- 
venient free parameter when using diagnostic diagrams to 
analyze emission-line ratios coming from different galaxies. 
However, in photoionization models for a given NLR where 
the distance D from the ionizing radiation source as well as 
the source luminosity are known, U is not a free parame- 
ter anymore but it is constrained by the source luminosity, 
the distance from the ionizing radiation source and the gas 
density (Ferguson et al. 1997). 

As shown in Tables 1, 2 and 3, the main emission line 
intensities were obtained as a function of the distance to 
the center. These measured line intensities result from the 
contribution of all the clouds in the line of sight covered by 
each aperture. It is not in the scope of this paper to provide 
a specific model for the observed data from each aperture. 
Our objective here is to give a scenario to try to explain the 
observations in general. Thus, in the following, we search for 
photoionization models which can reproduce the observed 
range of intensities. 

5.1 Multi-cloud photoionization model 

First attempts to model the emission-line region with a ho- 
mogeneous gas photoionized by the nuclear radiation failed 
to reproduce the observed NIR line intensities. As shown 
below, at least two density regions are necessary to explain 
the observed data. 

In the following, we present a two-cloud model for the 
NIR emission-lines using the photoionization code Aangaba 
(Gruenwald & Viegas 1992). The model must repro- 
duce the observed coronal lines [SIX]A1.252, [SVIII]A0.991, 
[SiX]A1.430, and [S VI]A1.962, as weU as the [SIII]AA0.907, 
0.953, and [SII]AA1.032 emission lines (Tables 1, 2 and 3). 

The sulfur line intensity ratios are not directly de- 
pendent on the adopted chemical abundance. Thus, they 
are first used to limit the other input parameters, i. e., 
U (ionization parameter), n (density) and D (distance 
from the ionizing source), given the ionizing radiation 
luminosity for NGC 1068, Lion- The model results are 
then compared to other line ratios, including those for 
[FeXI]A0.7889, [FeX]A0.6374, [FeVII]A0.6087 (Rodriguez- 
Ardila et al. 2006). We assume a solar chemical abundance 
(Grevesse & Anders 1989). 

For NGC 1068, the ionizing radiation luminosity is 
about lO'*'''^ erg/s assuming an UV power-law index equal 
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to -1.7 and an X-ray index of -0.5 (Pier et al. 1994, if the 
reflection factor is 0.01). 

In order to fit the high-ionization lines, a large value 
of U is usually required. However, due to the relationship 
between U, n, R and Lion, a high U requires a low n, and 
a low distance from the ionizing source. In the models pre- 
sented by the different authors referred above, the density 
value varies from 10^ cm^^ (Groves et al. 2006) to 10^ cm~"^ 
(Ozaki 2009), although some of them use filling factors as 
low as 6.5x10"*. 

From the observed coronal emission line intensities we 
note that for gas clouds extending to 200 pc, values of U > 
are required to fit high-ionization lines, and could be reached 
only for densities lower than 10^ cm~"^ (Rodriguez- Ardila et 
al. 2006, figure 12). After running a series of models for a 
homogeneous low density gas around the ionizing radiation 
source, we selected n = 1 cm"'' and U = 10 to better fit- 
ting the coronal lines; such a model is called model A. The 
line intensities relative to H/3 obtained for model A, for an 
optically thick cloud (r^^ — >■ oo), are given on the second 
column of Table 6. However, such a model also produces a 
large semi-ioruzed region, resulting in a [SIII]/[S II] line ra- 
tio too high. Thus the final model for the low density gas 
should be matter bound (see below). 

While the coronal line ratios are more sensitive to the 
U parameter, the [SIII]/[SII] line intensity ratio is more 
dependent on the gas density. The observed values indi- 
cate that densities 10^ cm""^ or higher are required to re- 
produced the observed line ratio. Notice that high density 
clouds (10^'^ cm~"^) are also favored by Ozaki (2009) when 
modeling the structure of the NLR of NGC 1068 using the 
optical emission-lines [FeVII], [OIII], and [OI], relative to 
H/3. Thus, we assume that high density clouds (10^ cm""') 
are radiation bound, located farther from the central source 
and being shielded from the UV, but not from the X-ray, by 
the low density gas producing the coronal lines. The results 
for a model using an ionizing spectrum with E ^ 200 eV are 
shown in Table 6 (model B). In addition to [SHI] and [SII] 
lines, this high density gas produces also the [SiVI] emission 
line, but not the other coronal lines as expected. 

The high density clouds contribute mainly to the low 
ionization emission lines. In order to obtain the fraction f of 
the contribution of the high density clouds relative to the 
coronal gas, we use the [SVIII]/[SIII] observed line ratios 
(Table 1, 2 and 3). Since [SHI] comes both from models A 
and B, such line ratio depends on the optical depth at the 
Lyman limit of the low density region. As noticed above, the 
[SiVI] line could partially come from the high density clouds. 
However their contribution is much smaller than that from 
the low density region, so that the [SiX]/[SiIV] ratio is prac- 
tically independent off. As said above, the low density cloud 
must be matter bound for not emitting too much [SII] or 
[SHI]. The best results are obtained with r^^ for model A 
between 1 and 5. In Table 6, the line intensity ratios for Tj^^ 
= 5 and 1 are listed (third and fourth columns, respectively). 
For T^y of model A equal 1 and 5, the corresponding frac- 
tions f are equal to 0.078 and 0.065, respectively, leading to 
a luminosity ratio of H/3 between the two components equal 
to 10.2 and 38.7. The range of the observed infrared line ra- 
tios is shown in Table 7, compared to the results of the two 
possible two-cloud models, i. e., the combination of model 
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Line intensities for components A and B, relative to 
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Ly 
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Ly 



mod A 

Ly 
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Ly 



[S n]/HI3 


3.5E-3 
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[S III] /H 13 


9.1E-1 
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[S Ylll]/Hl3 
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4.5E-2 
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[S IX]/HI3 
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3.5E-1 


O.OE+0 


[Si VI] /H/3 


6.6E-2 


9.1E-2 


2.5E-1 


4.2E-5 


[Si X]/ H/3 


2.1E-2 


2.9E-2 
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[Fe VII]////3 


l.lE-1 


1.5E-1 
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O.OE+0 


[Fe X]/Hfi 


4.9E-2 
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9.7E-2 
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Table 7. Result; 
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Ly 
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Ly 



[S III] /SI II] 


+0.8 to +1.2 


+1.0 


+1.1 


[S VIII] /[S III] 


-2.0 to -1.5 


-1.7 


-1.7 


[S IX] /[S VIII] 


+0.0 to +0.7 


+0.4 


+0.4 


[Si X]/[Si VI] 


-0.4 to +0.0 


-0.4 


-0.5 


[Fe X]/[Fe VII] 


-0.5 to -0 4 


-0.3 


-0.3 


[Fe XI] /[Fe X] 


-0.9 to -0.2 


-0.4 


-0.4 



B with model A (for r^^ =1 and 5). Notice that detailed 
models for each aperture may require a range of densities. 

All line ratios are in the range of the observational data 
(Tables 1, 2 and 3), including the coronal iron lines from 
Rodriguez- Ardila et al. (2006). In that paper, only models 
with density ^ 10^ cm""' were used, and pure photoioniza- 
tion models could not reproduce the line ratios. Here we 
relax this lower limit, and most of the ratios could be re- 
produced. This, however, does not discard the possibility of 
having other ionization mechanisms playing a role. To put 
a stronger constrain on this issue it would be necessary to 
fit a large amount of emission-lines in a broader range of 
wavelengths, which is out of the scope of this paper. 



6 THE CENTRAL MASS 

We have used the stellar absorption features to estimate the 
mass of the central region inside the extraction apertures. 
The velocity dispersion was estimated using the same veloc- 
ity distribution and methodology proposed by Van der Marel 
& Franx (1993). The h3 and /i4 coefficients, which measure 
the distortions from a gaussian distribution, are very uncer- 
tain and therefore we restrict our discussion to the velocity 
dispersions. The results were obtained using as template a 
mixture of standard stellar spectra with the same resolution 
as our nuclear spectra of NGC1068 and high S/N. The spec- 
tral region was the same used by Tamura et al. (1991) for 
this same galaxy, around 2.26 ]xm, where there are strong 
molecular absorptions due to CO. This region seems ade- 
quate since the interstellar absorption is small and there are 
no emission lines. Errors in the velocity dispersions were ob- 
tained using the S/N of the spectra. An example of the fit 
is shown in Figure 13. 
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Figure 14. Core model applied to NGC 1068 



Figure 13. Example of the fitting process. Color lines represent 
the observed spectra and dashed black lines are the ajusted model. 



The results are shown in Figure 14 and we can observe 
that the velocity dispersion rises strongly in the central re- 
gion, probably due to the AGN contribution. Moreover in 
this very central region the hght contribution of the AGN 
and the hot dust contaminates the stellar continuum and 
makes it difficult to measure the velocity dispersion. Away 
from the center the nuclear contamination is much lower 
in our extracted spectra and we observe that the velocity 
dispersion remains quite large for a typical spiral galaxy. 
Nevertheless moving outwards to 100 pc from the nucleus 
the velocity dispersion is quite close to the values already 
reported in the literature. 

Using these determinations and the Virial theorem we 
give a very crude estimate showing how the central mass de- 
pends on the adopted nuclear distance. Our approach con- 
sists in using a simple model that can be resolved analyti- 
cally and compared to the data. Basically the model consists 
of an stellar core with uniform density, p, , extending to a 
maximum radial distance R* . The total mass of this stellar 
core is M* and in the central region we add a punctual mass. 
Mo, that simulates the presence of the AGN and its accre- 
tion disk. The total mass profile can be readily estimated 
as 



M(7?) = M. (-|-) -t- 



Mn 



(3) 



where R represents the radial distance. Assuming that the 
velocity distribution of the stellar component is random and 
isotropic we derive that 



av{R)^GM, 



IR^ 



7 1 



5 i?3 + 3 ij 



(4) 



where (t^(R) is the velocity dispersion and 7 = Mo/M* is a 
measure of the relative importance of the central mass com- 
pared with the core stellar distribution. However we only 
have access to the average kinetic energy of the stars pro- 
jected in the line of sight of the observer. Weighting this 
component by its luminosity contribution, or by mass for a 
constant mass-luminosity ratio, we have 



2 ,, J^^ aliR)p{R)dl 
o- v.obs(r) = 



Jo PiR)dl 



(5) 



and assuming that the density of the stellar core is uniform, 
we finally obtain our estimation for the observed line of sight 
velocity dispersion as 



o- v,obs{r) 



^ V .ohs\-^* 



3 + 57 



(6) 



Ri ' {1-r^/Rl)^/^ \ r2 



1 + "^-^ + 



57 



1/2 



For large values of 7 the central mass dominates and we 
simply have a keplerian profile which is shown as the blue 
continuous curve in Figure 14. We can see that the observed 
profile is completely inconsistent using the usual BH mass 
estimation (~ 1 — 2xlO'^M0) found in the literature (Lodato 
& Bertin 2002 ). To explain the high velocity dispersion of 
255 km/s in the central region we need to include an stellar 
core of 1.5x lO^'^M© extending to approximately 200 pc from 
the nucleus shown as the dotted line. Therefore a value 7 = 
6.6 X 10"'* would be a more appropriated contribution of the 
central mass. 



T SUMMARY AND CONCLUSIONS 

We carried out long-slit spectroscopy in the Near-IR of 
NGC 1068 and presented here for the first time the whole 
range from 0.9 to 2.4[j.m, extending 15" in the nuclear re- 
gion. The nuclear emission shows many important emission 
lines, from low ionization to coronal lines, molecular lines 
and absorption features. We show that many emission lines 
are composed by two narrow components, one in the galaxy 
rest frame, and the other blueshifted in the south and red- 
shifted in the north apertures. The relative intensity of the 
peaks also varies and they seem to be correlated with the jet 
of the galaxy. No broad component (~ 2000) was identified 
in any line profile. Not all emission lines have double peaked 
profiles and some lines are present only near the center, par- 
ticularly the high ionization lines. This is a clear indication 
that the emitting gas is more ionized in the nuclear region. 
We estimated the extinction based on Paschen hydrogen 
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and [Fell] emission lines. E(B-V) is highest in the nucleus 
then in the adjacent regions, but there are peaks around 
150 - 200 pc from each side of the nucleus. This is also close 
to the region where the second components of the emission 
lines have a maximum, and may be associated with dense 
clouds around these regions. H2 emission is very extended, 
but it is difficult to identify the ionization process because 
only the l-O(Sl) line could be measured. 

Simple photoionization models were used to try to ex- 
plain the emission lines, and a first approximation shows 
that the line ratios can be explained with pure photoion- 
ization, as long as low density clouds are present. A more 
detailed analysis however is necessary, to rule out other ion- 
ization mechanisms. 

We also used the stellar signatures to determine the 
mass of the central region, obtaining a stellar core of 
1.5xlO^°M0. 
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